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Abstract: In arid desert regions of northwestern China, reclamation and subsequent irrigated cultivation 
have become effective ways to prevent desertification, expand arable croplands, and develop sustainable 
agricultural production. Improvement in soil texture and fertility is crucial to high soil quality and stable 
crop yield. However, knowledge on the long-term effects of the conversion of desert lands into arable 
croplands is very limited. To address this problem, we conducted this study in an arid desert region of 
northwestern China to understand the changes in soil physical-chemical properties after 0, 2, 5, 10, 17, and 
24 years of cultivation. Our results showed that silt and clay contents at the 17-year-old sites increased 17.5 
and 152.3 folds, respectively, compared with that at the 0-year-old sites. The soil aggregate size fraction and 
its stability exhibited an exponential growth trend with increasing cultivation ages, but no significant change 
was found for the proportion of soil macroaggregates (>5.00 mm) during the 17 years of cultivation. The 
soil organic carbon (SOC) content at the 24-year-old sites was 6.86 g/kg and increased 8.8 folds compared 
with that at the 0-year-old sites. The total (or available) nitrogen, phosphorus, and potassium contents 
showed significant increasing trends and reached higher values after 17 (or 24) years of cultivation. Changes 
in soil physical-chemical properties successively experienced slow, rapid, and stable development stages, but 
some key properties (such as soil aggregate stability and SOC) were still too low to meet the sustainable 
agricultural production. The results of this long-term study indicated that reasonable agricultural 
management, such as expanding no-tillage land area, returning straw to the fields, applying organic fertilizer, 
reducing chemical fertilizer application, and carrying out soil testing for formula fertilization, is urgently 
needed in arid desert regions. 
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1 Introduction 


Sustainable agricultural production and food security are important prerequisites to address 
challenges such as population growth, climate change, and land degradation (Lal, 2009; FAO, 
2017), particularly in arid desert regions that experience shortage of water resources, high 
temperature, and extreme and frequent droughts. In arid desert regions of northwestern China, 
agricultural reclamation through the conversion of desert lands to arable croplands has a long 
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history to resolve the conflicts between land desertification and agricultural development (Li et al., 
2006; Su et al., 2010). However, in recent decades, in addition to the expansion of oases, population 
growth, and extensive management of agricultural production, a number of environmental 
problems have emerged, including the drawdown of groundwater, salinization, misuse of chemical 
fertilizer, and abandonment of arable croplands, which limit sustainable agricultural development 
in these regions (Zhang et al., 2013). The protection and sustainable utilization of soil resources 
are the premise of sustainable agricultural production in arid desert regions. 

Soil texture and fertility in arid desert regions are characterized by loose structure, rich sand 
particles, weak biological effects, and low organic matter and nutrient levels (Verheye, 2009). 
However, soil texture and fertility can be significantly improved after reclamation due to irrigation, 
fertilization, cultivation, and other significant changes in the soil formation process, especially the 
increases of biomass and root carbon input and other organic matter input (Tiessen et al., 1992; 
Jaiyeoba, 2003; Caravaca et al., 2004; Su et al., 2010; Gelaw et al., 2014; El-Shahway et al., 2016; 
Xiao et al., 2019). Generally, the proportion of silt and clay contents in cultivated soils increases 
after reclamation, which is caused by the irrigation with silt-laden river water (Fullen et al., 1995), 
the increase in dust fall (Reynolds et al., 2001; Xu et al., 2016), and the decline in wind erosion 
(Duniway et al., 2019). Coupled with soil silt and clay particles, the decomposition of crop residues 
and organic fertilizer could help soil aggregate formation and soil organic carbon (SOC) 
sequestration in soils (Elliott, 1986; Lal, 2004; Wei et al., 2014). SOC in cultivated soils plays a 
vital role in improving soil fertility and quality via nutrient cycling, soil aggregation, wind erosion 
resistance, and agricultural production (Gelaw et al., 2014; Xiao et al., 2020). Moreover, SOC 
represents a large carbon pool despite its low contents. A previous study estimated that the world's 
dryland soils contain 241 Pg of SOC, which was approximately 40 times the CO2 emissions caused 
by human activities during the 1990s (Lal, 2004). Agricultural activities that convert desert lands 
into arable croplands likely accelerate the carbon sequestration of arid ecosystems. 

The recovery of soil texture and fertility is a long-term and complex process that depends on the 
soil parent materials, climate, tillage methods, and agricultural management. Xiao et al. (2009) 
showed that 10 years of cultivation significantly increased the SOC pool but not the 
macroaggregation of desert soils. Su et al. (2010) found that absolute amounts of soil aggregates, 
SOC, and nutrients were too low to support sustained crop yield after 40 years of cultivation. 
However, existing studies are insufficient to fully understand changes in soil texture and fertility 
after continued cultivation, especially in arid desert regions of northwestern China. The objective 
of this study is to test how reclamation and subsequent irrigated cultivation influence soil texture 
and fertility during a long-term cultivation period in an arid desert environment in northwestern 
China. To achieve this goal, we systematically assessed changes in soil physical-chemical 
properties in response to cultivation ages, the contribution of principal components of soil 
properties, and the development tendency of soil succession. This study is helpful for rational 
reclamation and subsequent irrigated cultivation of desert lands and scientific management of soil 
fertility at the expanded arable croplands. 


2 Materials and methods 


2.1 Study area 


The Aral reclamation region (40°20'16"—41°47'18"N, 79°22'33"—81°53'45"E; 996-1097 m a.s.l.), 
located in Xinjiang Uygur Autonomous Region of northwestern China, is a representative arid 
desert region with extreme drought in the world (Fig. 1). As anew center of agricultural production 
since 1958, the Aral reclamation region occupies the alluvial plain of the Aksu River and Yarkant 
River in the upper reaches of the Tarim River. This region has a typical desert climate and is 
characterized by dry hot summers and cold winters, with an annual mean temperature of 10°C. The 
mean annual precipitation is 67.2 mm, which is mainly concentrated from June to August, 
accounting for more than 60.00% of the annual precipitation. The mean annual potential 
evaporation reaches 2110.5 mm and the annual crop evaporation during the growing period from 
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April to October is 1797.0 mm, occupying 85.19% of the annual evaporation. The mean annual 
sunshine hours are 2950 h and the frost-free period is 207 d. The zonal soil is meadow solonchak 
derived from diluvial-alluvial materials of the Tarim River. Because of population growth and the 
need for agricultural development, desert lands have been gradually converted into arable croplands 
with irrigation since the 1950s. Therefore, this region becomes the largest cotton (Gossypium 
hirsutum L..) production base and the only long-staple cotton production base in China. 
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Fig. 1 Geographical location of the study area in Xinjiang Uygur Autonomous Region of northwestern China (a) 
and overview of the sampling locations (b) from the Google Earth images. Numbers ®-@© represent arable 
croplands with 0, 2, 5, 10, 17, and 24 years of cultivation, respectively. 


2.2 Experimental design and soil sampling 


In this study, we adopted the space for time substitution method to choose croplands with 
different cultivation ages (including 2, 5, 10, 17, and 24 years; Fig. 1). The adjacent natural desert 
land (0 year of cultivation) was used as a control for the comparison with the cultivated croplands 
(2, 5, 10, 17, and 24 years of cultivation). The major crop grown in the above cultivated croplands 
is cotton, with only one harvest annually. The annual fertilization rates of cotton are 
approximately 360-460 kg N/hm?, 90-140 kg P20s/hm?, and 50-90 kg K2O/hm?. The total 
irrigation amount is approximately 7300-10,200 m°/hm? with flooding irrigation in the off season 
and drip irrigation in the growing season. 

In mid-September 2017, six replicate sites were randomly selected within each cropland 
sampling location (including the control). At each site, soil samples from soil depths of 0-50 cm 
were collected by a soil auger and mixed as a composite sample. Overall, 36 composite soil samples 
were collected (six cropland sampling locationsxsix replicate sites) and transported to the 
laboratory. After air drying, one aliquot of each composite soil sample was passed through a 2.00- 
mm sieve to analyze the soil particle size distribution. These soil samples were also further ground 
and sieved at 0.25 mm to analyze the following properties: SOC, total nitrogen (TN), total 
phosphorus (TP), total potassium (TK), available nitrogen (AN), available phosphorus (AP), and 
available potassium (AK). The other aliquot of each composite soil sample was passed through a 
5.00 mm-sieve to analyze the soil aggregate separation. 


2.3 Measurements of basic soil properties 


Soil bulk density was measured with the cutting ring method in situ. The particle size fraction of soil 
samples was determined by a laser diffraction particle size analyzer (S3500, Microtrac Inc., USA). 
Soil pH was measured with a soil to water ratio of 1:2.5, and the total salt content was determined by 
the residue drying method. SOC was measured with the Walkley and Black's dichromate oxidation 
method (Nelson and Sommers, 1982). TN was determined via the micro Kjeldahl digestion procedure 
(Foss Kjeltec 2300, Denmark) and AN was measured by the alkaline diffusion method. TP was 
analyzed colorimetrically by a spectrophotometer after soil digestion with HC]O4,-H2SO, and AP was 
determined by the Olsen method. TK and AK were measured by the flame photometry after HF— 
HClO; digestion and the colorimetric method after NH4OAc extraction, respectively. 


2.4 Determination of soil aggregate size fraction and its stability 
Soil aggregate size was measured by the wet sieving apparatus (Li et al., 2006). A total of 50 g of 
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air-dried soil samples was transferred to a nest of five sieves with square openings of 0.25, 0.50, 
1.00, 2.00, and 5.00 mm. Then, the sieves were put into a bucket and deionized water was poured 
along the wall of the bucket until the soil aggregates submerged in the water. After 5 min of 
moistening and soaking, the bucket was oscillated for 5 min at an amplitude of 3 cm and a frequency 
of 30 cycles/min. The remaining soil sample on each sieve was dried at 60°C for 48 h to a constant 
weight. Soil aggregate size distribution was expressed by the mean weight diameter (MWD; mm) 
of aggregate (Elliott, 1986): 


6 
MWD = >" DW, (1) 
i=1 


where D; is the average diameter of soil aggregate between adjacent size classes (mm); and W; is 
the mass percentage of the remaining aggregate on the i" sieve (%). 


Water stable aggregate percentage (WSAP; %) was calculated as follows: 


WSAP = | es ] x100%, (2) 


where Mwsa is the corrected mass of water stable aggregate with diameter greater than 0.25 mm (g); 
and M. is the adjusted dry mass of soil aggregate with diameter of 2.00-5.00 mm (g). 


2.5 Statistical analysis 


In this study, the descriptive statistics (mean and standard error) for soil physical-chemical 
properties in croplands with different cultivation ages were calculated. One-way ANOVA was 
performed to test the effects of cultivation age on soil physical-chemical properties. When the 
above effects were significant (P<0.05), Tukey's HSD post hoc test was performed to compare the 
mean values of above soil physical-chemical properties in croplands with different cultivation ages. 
Regression analyses were used to indicate the trends of soil physical-chemical properties along the 
cultivation time sequence. Principal component analysis (PCA) was conducted to characterize the 
associations of soil physical-chemical properties in croplands with different cultivation ages. 
Hierarchical clustering was applied to the soil physical-chemical property data using the Ward's 
method to indicate the different recovery stages of soil environment. All of the statistical analyses 
were conducted using JMP software v.10.0.0 (SAS Institute, Cary, NC, USA). 


3 Results 


3.1 Changes in soil particle size fraction and bulk density 


The soils in the study area were originated from both the alluvial soils and deposited aeolian sand, 
which were mainly composed of medium-fine sand with the content ranging between 95.95% and 
96.85% (Table 1). After reclamation and subsequent agricultural activities, the sand content 
decreased from 96.40% to 31.60% for the 17-year-old sites, while the contents of silt and clay 
increased 17.5 and 152.3 folds compared with those at the 0-year-old sites (control; Table 1), 
respectively. Interestingly, we found that soil particle size fraction reached the extreme values at 
the 17-year-old sites. Specifically, the sand content had a minimum value of 31.60%, and the silt 
and clay contents had the maximum values of 62.27% and 6.09%, respectively. Moreover, a 
significant change in particle size fraction only occurred in soils with more than 5 years of 
cultivation. Soil bulk density generally fluctuated as the increasing of cultivation ages. After 24 
years of cultivation, soil bulk density reached the minimum value of 1.44 g/cm? and decreased by 
4.55% compared with that at the 0-year-old sites. 


3.2 Soil aggregate size fraction and its stability 


As cultivation prolonged, soil aggregate size fraction changed significantly among croplands 
with different cultivation ages. Generally, the proportion of soil aggregate sizes (i.e., 0.25—0.50, 
0.50—1.00, 1.00-2.00, 2.00-5.00 and >5.00 mm) demonstrated an increasing trend and significantly 
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reached the highest values at the 24-year-old sites (Table 1). Due to the limitation of original soil 
conditions with the shortage of binding agents such as organic matter and clay, the proportion of 
large macroaggregate (>5.00 mm) had no significant differences among the 2-, 5-, 10-, and 17- 
year-old sites. After 24 years of cultivation, soil aggregate size fraction mainly contained the 
microaggregate (0.25—0.50 mm) with the proportion of 15.06%. MWD showed significant changes 
with cultivation ages, ranging from 0.04 mm at the 2-year-old sites to 0.39 mm at the 24-year-old 
sites (Table 1). Changes in MWD exhibited an exponential growth trend along the cultivation 
chronosequence (Fig. 2). Similarly, WSAP increased significantly with increasing cultivation ages 
and exhibited an exponential growth trend (Fig. 2), i.e., increasing from 7.42% at the 0-year-old 
sites to 42.40% at the 24-year-old sites (Table 1). 
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Fig. 2 Exponential relationships between sites with different cultivation ages and MWD (a) and between sites 
with different cultivation ages and WSAP (b). MWD, mean weight diameter; WSAP, water stable aggregate 
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3.3 Soil chemical properties 


After reclamation and subsequent cultivation activities, soil chemical properties were ameliorated 
gradually with increasing cultivation ages. SOC content at the 24-year-old sites was significantly 
higher than those at the other sites, and it increased 8.8 folds more than that at the 0-year-old sites 
(Table 1). Contents of TN and TP at the 17- and 24-year-old sites reached higher levels compared 
with those at the other younger sites. AN exhibited a unimodal pattern with the cultivation age 
sequence and reached its peak value at the 17-year-old sites. AP showed a significant increase after 
24 years of cultivation. AK reached its highest value at the 17-year-old sites, i.e., 3.3 times higher 
than that at the 0-year-old sites. The C/N ratio showed a significant increase trend with increasing 
cultivation ages and the values at the 2-, 5-, 10-, 17-, and 24-year-old sites were 5.5, 2.7, 2.6, 3.7, 
and 2.9 times higher than that at the 0-year-old sites, respectively. A significant decline in soil pH 
was detected after 2 years of cultivation. Moreover, after 24 years of cultivation, soil pH decreased 
by 7.1% in comparison with that at the 0-year-old sites (Table 1). 


3.4 Multivariate statistical analysis results 


PCA identified three components for all the cultivated and uncultivated sites, and these three 
components explained 81.6% of the total variation. The first component (PC1) explained 40.6% of 
the total variation, and its loadings were soil bulk density, soil aggregates of different size fractions 
(0.25, 0.50, 1.00, 2.00, and 5.00 mm), MWD, WSAP, SOC, TN, TP, and AP (Table 2; Figs. 3a and 
b). The second component (PC2) explained 32.6% of the total variation and was loaded by clay, 
silt, sand, TK, AN, and AK (Table 2; Figs. 3a and b). Soil pH was loaded on the third component 
(PC3), which explained 8.4% of the total variation (Table 2; Figs. 3a and b). Across the PC1 axis, 
the cultivated and uncultivated sites were arranged in order with cultivation ages, namely, the 0-, 
2-, 5-, and 10-year-old sites on the negative direction of the PC1 axis which can be divided into a 
group (Fig. 3c), and the 17- and 24-year-old sites on the positive direction of the PC1 axis with a 
decentralized distribution (Figs. 3c and d). 
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Table2 Correlation coefficients between soil physical-chemical properties and the three principal component 
analysis (PCA) axes and the total variation explained by the three PCA axes 


Soil physical-chemical 


ae TEES PC1 PC2 PC3 
Clay 0.301 0.925 —0.146 
Silt 0.363 0.876 —0.181 
Sand —0.358 —0.882 0.178 
Bulk density —0.496 0.095 —0.202 
A0.25 0.871 0.316 —0.043 
A0.50 0.875 0.062 0.033 
A1.00 0.848 0.173 0.313 
A2.00 0.889 0.136 —0.187 
A5.00 0.598 0.442 0.366 
MWD 0.942 0.194 —0.133 
WSAP 0.955 0.200 —0.102 
SOC 0.773 0.442 —0.188 
TN 0.617 0.524 —0.069 
TP 0.621 0.559 0.292 
TK 0.023 0.904 0.121 
AN 0.372 0.799 —0.237 
AP 0.861 0.048 —0.063 
AK —0.080 0.945 —0.037 
pH —0.265 —0.207 0.838 
Total variation explained 40.6% 32.6% 8.4% 


Note: A0.25, A0.50, A1.00, A2.00, and A5.00 refer to soil aggregate size fractions of 0.25-0.50, 0.50-1.00, 1.00-2.00, 2.00-5.00, 
and >5.00 mm, respectively. PC1, PC2, and PC3 refer to the first, second, and third components, respectively. 


Hierarchical cluster analysis with the Ward's method was employed to identify similarity of soil 
properties for the cultivated sites, which rendered a dendrogram that grouped the cultivated and 
uncultivated sites into three statistically significant clusters when the Ward's distance was 7.62 (Fig. 
4). Cluster 1 included the 0-, 2-, 5-, and 10-year-old sites that represented the 1* reclamation stage 
with slow development of soil properties; cluster 2 (17-year-old sites) corresponded to the rapid 
amelioration of soil conditions, which was classified as the 2"d reclamation stage; and cluster 3 
included the 24-year-old sites and represented the stable soil improvement stage (Fig. 4). The three 
clusters corresponded to the biplot of PCA (Fig. 3c). 


4 Discussion 


4.1 Impact of reclamation on soil physical properties 


Continued cultivation is an effective way to improve the particle size fraction of arid desert soils. 
Our results found that silt and clay contents increased significantly while sand content decreased 
after reclamation and subsequent irrigated cultivation (Table 1). Originally, the soil parent materials 
in our study area come from the weathering materials of the parent rock of the Tianshan Mountains 
in quaternary sediments, which are deposited by the Tarim River and Yarkant River. At the same 
time, under the influence of the Taklamakan Desert, the aeolian sediments occur in the reclamation 
region, and aeolian sediments and alluvial sediments often alternately participate in soil formation 
(Gong et al., 2016). Jaiyeoba (2003) observed that cultivated sites had higher silt content than 
barren land. A previous study conducted in an arid oasis showed that silt and clay contents increased 
from 4.50% to 20.90% and from 4.80% to 8.90%, respectively, after approximately 40 years of 
continued cultivation (Su et al., 2010). Due to limited precipitation, irrigation water from the inland 
Tarim River, Aksu River and other reservoirs ensures the sustainable development of agriculture in 
the study area. According to the official statistics, the Tarim River has an annual average sediment 
(mostly silt and clay) concentration of 4.65 kg/m? and a maximum sediment concentration of 20.9 
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Fig. 3 Principal component analysis (PCA) showing loading values of soil physical-chemical properties for the 
first component (PC1) and the second component (PC2) axes (a) and the PC1 and the third component (PC3) axes 
(b); and the score plots of the cultivated and uncultivated sites with different cultivation ages along with the PC1 
and PC2 axes (c) and the PC1 and PC3 axes (d). Error bars show the standard errors (n=6). BD, bulk density; SOC, 
soil organic carbon; TN, total nitrogen; TP, total phosphorus; TK, total potassium; AN, available nitrogen; AP, 
available phosphorus; AK, available potassium. A0.25, A0.50, A1.00, A2.00, and A5.00 refer to soil aggregate size 
fractions of 0.25-0.50, 0.50-1.00, 1.00-2.00, 2.00-5.00, and >5.00 mm, respectively. 


kg/m? (Cai, 2013). All these silt-laden river water contributes to the increase in silt and clay 
contents, similar to some oasis croplands in arid desert regions of northwestern China (Fullen et 
al., 1995; Lal, 2004). Under the protection of shelterbelts and crops, the wind erosion of surface 
soils decreased significantly while the atmospheric dust increased. For example, Fearnehough et 
al. (1998) reported that the annual dustfall in arid desert regions ranges from 21.5 to 56.0 mg/cm?, 
which promotes the silt and clay contents in the surface soil layer of the cultivated regions. 

The soil aggregate size fraction showed a similar tendency as the soil particle size fraction, 
while the regression analysis demonstrated a positive linear relationship of WMD and WASP with 
silt and clay contents, and a negative linear relationship among WMD, WASP, and sand content 
(Fig. S1). Therefore, increases in silt and clay contents can be considered a prerequisite of soil 
aggregate formation (Li et al., 2006; Verheye, 2009; Su et al., 2010). In addition to soil minerals 
(clay particles, silt, and fine sand), plant roots, fungal hyphae, remains, bacteria, and free 
amorphous organic matter are also important contributors to soil aggregates (Tisdall and Oades, 
1982; Tiessen et al., 1992; Six et al., 2000; Ouyang et al., 2018). Due to the extremely weak 
biological effect on soil formation, the uncultivated soils have very low organic matter and nutrient 
contents., Irrigation and agricultural fertilizer application after reclamation increase the input of 
organic materials and continuous agricultural utilization improve the input of biomass to the soils, 
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Fig. 4 Dendrogram results from the Ward's method of hierarchical cluster analysis for the cultivated and 
uncultivated sites. The values of 5.05, 7.62, and 10.27 refer to the rescaled distances for the different cluster 
combinations. Each cropland sampling location (including the control) had six replicates. The curve with crosses 
in the figure represents the stages of cluster combined. 


therefore accelerating the formation of soil aggregate (Tiessen et al., 1992; Hannachi et al., 2015). 
In addition, microbial partners, such as Arbuscular mycorrhizae fungi, associated with crop roots 
synthesize mucigels and polysaccharides that bind and transform the microaggregates into steady 
macroaggregates (Tisdall et al., 1997). Our results indicated that after 24 years of cultivation, soil 
aggregate size fraction contained the microaggregate (0.25-0.50 mm) with the proportion of 
15.06% and the macroaggregate in a low level with the proportion of 0.50% (Table 1). One reason 
for this scenario is that low SOC concentration is not beneficial for binding isolated sand particles 
together with very limited silt and clay contents (Li et al., 2006; Su et al., 2010) (Fig. S2). Another 
reason is that macroaggregates are less stable than microaggregates and sensitive to agricultural 
practices, e.g., plowing brings soils from the deep layer to the surface layer and physically breaks 
down macroaggregates, resulting in an increas in the turnover of SOC (Caravaca et al., 2004; Wei 
et al., 2014). 


4.2 Impact of reclamation on soil chemical properties 


Our study showed that cultivation of desert lands accelerated the increase in soil nutrient properties, 
the tendency of which was consistent with the changes in soil particle size and aggregate size 
fractions (Table 1). A number of previous studies have indicated that agricultural practices (e.g., 
irrigation, tillage, fertilization, and crop rotation) could improve soil carbon sequestration in arid 
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and semi-arid regions after conversion of desert lands into arable croplands (Lal, 2004; Li et al., 
2006; Wu et al., 2008; Li et al., 2010; Hannachi et al., 2015). The potential increase in SOC and 
nutrient contents mostly depends on the initial fertility level of the uncultivated land. Su et al. 
(2010) found that soils with a low initial SOC level have a greater potential for carbon sequestration 
and TN increase. Based on a long-term study of irrigated agricultural soils in California, Wu et al. 
(2008) revealed that after 85 years of cultivation, irrigated farming significantly increased SOC and 
soil inorganic carbon through the addition of crop residues, input of root mass, and amendment of 
organic fertility (Dong et al., 2018; Medeiros et al., 2020) or through the conversion of soil CO2 to 
carbonates (Entry et al., 2004). Irrigated cultivation in arid and semi-arid desert regions may 
provide a path to contribute to carbon sequestration and reduce atmospheric CO2 concentrations. 
For example, Lueking and Schepers (1985) estimated that irrigated cultivation could increase SOC 
by 1.66 Mg/hm? within the top 30 cm of soil after 15 years of cultivation. In addition, due to 
microbial SOC decomposition, irrigated cultivation can also change the magnitude of soil CO2 
efflux and the input and output of soil carbon pool, and even influence the soil carbon balance of 
cultivated croplands (Li et al., 2010). The current study indicated that after 24 years of cultivation, 
SOC continued to rise and the magnitude of soil carbon sequestration was higher than that of carbon 
loss. 

Other soil chemical properties, including TN, TP, and TK, showed a similar trajectory as SOC 
during the different cultivation periods, and SOC had a positive correlation with TN, TP, and TK 
(Fig. S2). Under irrigated conditions, nitrogen (N) and phosphorus (P) are important factors 
restricting agricultural productivity and mostly rely on the application of chemical fertilizers. 
According to our field investigation, the actual application rates of N, P, and potassium (K) were 
264, 69, and 82 kg/hm?, respectively, which exceeded the recommended fertilization standards 
(230, 60, and 72 kg/hm?, respectively) for cotton cultivation in Xinjiang Uygur Autonomous 
Region, China. In fact, intensive high-yield agriculture is dependent on the addition of fertilizer. 
From 1960 to 1995, the global use of N and P increased 7.0 and 3.5 folds, respectively, and both 
are expected to increase another 3.0 folds by 2050 unless there is an improvement in fertilizer 
efficiency (Tilman et al., 2001). Undoubtedly, the conversion of desert lands into arable croplands 
contributed to the improvement in fertilizer efficiency because we found that increasing soil N, P, 
and K were positively correlated to soil physical proporties (including silt and clay contents, and 
MDW) in sites with different cultivation ages (Fig. S2). However, excess application of N fertilizer 
resulted in a low C:N ratio of 8.16 after 24 years of cultivation (Table 1), which was lower than the 
average values of soils in China and in the world (11.38 and 13.33, respectively; Aitkenhead and 
McDowell, 2000). 


4.3 Soil fertility and recovery potential 


Soil physical-chemical properties are the key factors in determining the soil fertility. According to 
the results of PCA, the PC1 was related to the "soil aggregate stability and chemical properties" set 
and was loaded by soil bulk density, soil aggregate size fractions, MWD, WSAP, SOC, TN, TP, and 
AP (Table 2). Soil aggregate stability, as an indicator of soil structure, resulted from the 
rearrangement of particles, flocculation and cementation (Duiker et al. 2003) and was directly 
related and mediated by SOC, biota, ionic bridging, and other soil nutrients (Fig. S2). Previous 
studies have shown that the water stability of soil aggregates depends on soil organic matter, which 
acts as a persistent organic binding agent and as a core in the formation of aggregates (Tisdall and 
Oades, 1982; Bronick and Lal, 2005). In the arid desert region of this study, due to irrigation, 
fertilization, cultivation, and other significant changes in the soil formation process, MWD, WASP, 
and soil chemical properties all showed an exponential growth trend during the 24 years of 
cultivation (Fig. 2). The PCA results indicated that soil particle size fraction, TK, AK, and AN were 
loaded on the second PCA axis, which was mainly associated with soil texture and available 
nutrients. During the cultivation process, increasing silt and clay contents play a vital role in the 
formation of favorable structures of sandy soils (Mitchell et al., 1998; Su et al., 2010). In addition, 
most soil K is oriented from K-bearing primary materials, such as muscovite, feldspars, and 
muscovite (Zörb et al., 2014). Both water-soluble and exchangeable K which are easily available 
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for crops, are related to clay content in soils (Barré et al., 2008). Our results also indicated that clay 
and silt contents were positively related to TK (Fig. S2). The third PCA axis was loaded by soil pH, 
which decreased by 7.1% after 24 years of cultivation in comparison with that at the 0-year-old 
sites (Table 2; Figs. 3a and b). In arid desert regions, soil salinization is one of the important factors 
restricting agricultural production, so a decrease in soil pH is beneficial to crop growth. High levels 
of N fertilization in our study area could directly or indirectly contribute to the decline of soil pH 
(Guo et al., 2010). 

Our results indicated that soil physical-chemical properties in sites with different cultivation ages 
could be divided into three stages, i.e., the slow growth period of the 1* reclamation stage, the rapid 
amelioration period of the 2™ reclamation stage, and the stable improvement period of the 3'4 
reclamation stage (Fig. 4). After acclimation, under the effects of irrigation, fertilization, and 
cultivation, the input of crop root biomass and other organic materials slowly improve the soil 
fertility conditions. In arid desert regions, due to the long-term wind erosion, desert soils are rich 
in sand particles, with loose structure and weak biological effects on soil formation; therefore, the 
process of improving soil texture and fertility undergoes long-term and complex periods. Study on 
the oasis arable croplands of the Hexi Corridor in northwestern China showed that the absolute 
amounts of soil aggregates, SOC, and nutrients were insufficient to maintain sustainable 
agricultural production after 40 years of cultivation (Su et al., 2010). Li et al. (2006) pointed out 
that aperiodic wind erosion conditions were not conducive to the stability of soil aggregates and 
the maintenance of soil fertility. Based on data from both the soil fertility investigation (2011— 
2012) and the second national soil survey (1980s), the soil fertility of cultivated land in the Aral 
reclamation region was relatively lower in organic matter and higher in alkali hydrolyzed nitrogen 
and AP (Cai, 2013). 


5 Conclusions 


From the perspective of soil physical-chemical properties, this study showed the development of 
soil texture and fertility and concluded that the conversion of desert lands into arable croplands 
under long-term cultivation improved the soil physical-chemical properties. Specifically, soil 
particle size fraction, soil aggregate stability, SOC, TN, AN, TP, AP, TK, and AK exhibited 
sustained growth trends with increasing cultivation ages. However, our data showed that the 
process of improving soil texture and fertility undergoes long-term and complex periods, because 
after 24 years of cultivation, some key soil physical-chemical properties (such as soil aggregate 
stability, SOC, etc.) were still too low to meet sustainable agricultural production. A systematic and 
comprehensive assessment, such as the soil quality index, is needed in the future to evaluate 
reclamation-induced changes in soil quality or fertility. For the future agricultural practices to 
effectively prevent wind erosion of arable croplands in arid and semi-arid desert regions, a range 
of agricultural management measures will be needed, such as expanding no-tillage land area, 
returning straw to fields, applying organic fertilizer, reducing chemical fertilizer, and carrying out 
soil testing for fertilization formulas. 
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Fig. S1 Linear relationships between soil aggregate stability (MWD and WSAP) and soil particle size fractions 
(sand, clay, and silt) of the cultivated sites. MWD, mean weight diameter; WSAP, water stable aggregate 
percentage. “* indicates significant correlation at P<0.001 level. 
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Fig. S2 Correlation coefficient matrix of soil physical-chemical properties. The red line denotes the linear 
relationship between each of two soil properties, and the red shadow denotes the confidence level of 95%. 


